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Abstract 

A feasibility study has been initiated on the use 
of remote earth observations for augmenting stream 
discharge prediction for the design and/or operation 
of major reservoir systems, pumping systems and 
irrigation systems. The near-term objectives are 
the interpolation of sparsely instrumented precipi- 
tation surveillance networks and the direct measure- 
ment of water loss by evaporation. The first steps 
of the study covered a survey of existing reservoir 
systems, stream discharge prediction methods, 
gage networks and the development of a self-adaptive 
variation of the Kentucky Watershed model, SNOPSET, 
that includes snowmelt. As a result of these studies, 
a special three channel scanner is being built for a 
small aircraft, which should provide snow, tempera- 
ture and water vapor maps for the spatial and tem- 
poral interpolation of stream gages. The reservoir 
system of the Western Division of the Bureau of 
Reclamation was chosen for future demonstration of 
how such remote observations might augment stream 
discharge estimates. 
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Introduction 

Remote earth observations from aircraft are 
being utilized to study the feasibility of applying 
future space payloads to stream discharge predic- 
tions. Such predictions would help in breaking the 
cycles between floods and droughts by distributing 
water more uniformly throughout the year. The 
more uniform distribution will preserve fertile 
lands and improve the production of food and fiber. 
The associated management of agricultural and 
forest resources follows directly from the water 
resources management. 

Damaging floods in the Mississippi and Tennes- 
see River basins were frequent before the develop- 
ment of water resource management systems by 
the Corps of Engineers and the Tennessee Valley 
Authority. Similar systems are now needed for 
other major river systems in underdeveloped coun- 
tries. The Department of Civil Engineering at 
Colorado State University has worked on such devel- 
opment projects for several years with the govern- 
ments of East Pakistan, West Pakistan, India, 
Thailand, and Venezuela. 

The following paper describes the first steps 
of a feasibility demonstration project. The long- 
range objective of this project is to integrate sta- 
tistical methods for machine interpretation of earth 
observations with hydrological simulation models 
for predicting the stream discharge of large drainage 
basins for the design and operation of major reser- 
voir systems, pumping systems, and irrigation 
systems. 

The short-term objective is to utilize multi- 
spectral observations of the Colorado River basin 
that will be obtained with MSFC instrumented air- 
craft for the following two applications: 
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1. The interpolation of sparsely instrumented 
precipitation surveillance networks with remote 
surveys 

2. The direct measurement of the water loss 
by evaporation. 

Hydrological Models 

Throughout the history of science the develop- 
ment of improved measurement techniques has 
stimulated the improvement of theoretical models 
of our physical environment. That is the situation 
with respect to remote sensing at this time. One 
very important area in which improved measure- 
ment techniques is stimulating the development of 
new models is that of water resources. Modern 
water resources management is based on hydrologi- 
cal simulation models. These fall into two categories: 
the first is statistical hydrology while the second 
consists of parametric hydrological simulation models. 
Statistical hydrological models may be used where 
many years of records are available but may not be 
readily extended to regions where historical records 
are not available. 

Parametric simulation models may be used in 
regions where recorded streamflow data are not 
available and may be used to assess the effects of 
changes in watersheds as well as to extend records 
in regions where streamflow or precipitation records 
are inadequate. 

Stanford University has developed one para- 
metric model which is sufficiently general for appli- 
cations to a great many different regions [ 1 ] . This 
model is illustrated in Figure 1. It may be described 
as a set of transfer functions which relate precipita- 
tion gage readings to stream gage readings. The 
precipitation gage readings are usually supplied from 
a measurement network within a drainage boundary. 
The standard format of the Bureau of Reclamation 
calls for hourly or daily readings of precipitation. 
These readings refer to the height of water or snow 
depth. During or shortly after rainfall or snowfall, 
the precipitation readings will give the direct input 
of available water. Measurements of evaporation 
from pans filled with water are used to indirectly 
estimate the water loss from wet soils and snow and 
by transpiration from plants. This loss is inferred 
by accounting for the changes of the measured evap- 
oration due to soil materials, surface slopes and 
vegetation cover. In most drainage basins such actual 
evaporation losses exceed the stream discharge. 


A successful hydrological simulation model 
should provide a reasonable agreement between the 
stream discharge that is predicted from the precipi- 
tation readings and the actual discharge readings of 
the stream gages. Such success has been achieved 
only for small and well instrumented drainage basins. 
The runoff prediction of larger watersheds is diffi- 
cult for two basic reasons. The first is the inade- 
quate coverage of larger watersheds with precipita- 
tion gages. A successful simulation model requires 
a precipitation surveillance network which covers 
various infiltration conditions (interflow) , slopes, 
and vegetation covers. The second reason is the 
uncertainty in the accurate prediction of evapotran- 
spiration losses from indirect measurements. 

Hydrological simulation models are needed to 
improve the operation of existing water resource 
management systems in large drainage basins and to 
design future reservoirs, irrigation systems, and 
cloud seeding operations for the more even distribu- 
tion of water throughout the year. The need to im- 
prove the prediction of stream discharge for existing 
systems in the U.S. may be illustrated for the inte- 
grated water supply and utilization system of the 
Western Division System of the U.S. Bureau of 
Reclamation, which operates 22 reservoirs, 16 
power plants, and 3 pumping stations (Fig. 2). Be- 
cause of an inadequate precipitation network in an 
isolated region, the predictions of the Sweet Water 
River discharge, which originates in the Wind River 
Range, were off by 800 percent in 1969. This pre- 
diction would have led to extensive flood damage in 
Casper, Wyoming, if the reservoirs above Casper 
had been filled close to capacity early in the snow- 
melt season. 

In mountainous regions, such as Colorado and 
Wyoming, the primary input for streamflow is winter 
precipitation and snowmelt. The detailed processes 
which produce streamflow from snowmelt are not 
well understood at this time [2 i . However, many 
past investigations [3 through 10] have indicated 
that every watershed has characteristic relationships 
between snowpack depletion and streamflow runoff. 
Relatively simple relationships can be derived and in- 
corporated into a model that may be used to synthe- 
size streamflow in any given basin. 

The Stanford Watershed Model [ 1 ] has been used 
for this purpose, but it is necessary to obtain the 
parameters through a trial-and-error adjustment 
procedure. This is not satisfactory since such a sub- 
jective approach makes extrapolation to undeveloped 
regions difficult. The model parameters obtained in 
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this manner will be different depending upon the in- 
dividual investigator and his understanding of the 
various phases of the hydrological sequence. This 
problem can be corrected by the use of a self- 
calibrating model which adjusts the parameters 
based upon a quantitative figure of merit. 

A self -calibrating model called OPSET (Opti- 
mal Set of Parameters) [ 11 ] has been developed 
from an extension [ 12 ] of the Stanford Watershed 
Model. However, this extension does not allow for 
snowmelt [ 13 ] , which is the primary input for 
streamflow runoff in the mountainous regions of 
interest. Another attempt at developing a self- 
calibrating watershed model was made at the Univer- 
sity of London in 1970 [ 14] . This was also based on 
the Stanford Watershed Model and a numerical optimi- 
zation technique developed by Rosenbrock [ 15 ] . The 
Rosenbrock technique was applied to the Kentucky 
Watershed Model by Colorado State University in- 
vestigators but it was found to be entirely too con- 
suming of computer time (Control Data Corporation 
6400 computer) and consequently a gradient proce- 
dure was developed and incorporated into the Ken- 
tucky Watershed Model. This self-calibrating model 
includes snowmelt and has been designated as 
SNOPSET. The first results of this model are given 
here (pp. 112-113) and provide the basis for the 
approach toward utilization of remote observations. 

The Existing Surveillance Network 

The feasibility of discharge prediction in un- 
developed regions might be demonstrated by expand- 
ing self-calibrating hydrological models from well 
to less well instrumented watersheds. The reservoir 
system of the Bureau of Reclamation provides such 
an opportunity as shown in Figure 2. Most of the 
region has only very few precipitation and stream- 
flow gages, which have, however, been recorded 
for many decades. These few gages support the 
present use of statistical hydrological models for 
the operational discharge forecast. However, the 
Bureau has now a pilot program of increasing the 
winter snowpack (Project Skywater) by cloud seeding 
by at least 30 percent. If initiated, such change of 
precipitation patterns might decrease the validity of 
the present statistical forecast, and parametric 
models of the SNOPSET type might thus be considered. 

A new precipitation and streamflow gage net- 
work is presently being installed by Western Scien- 
tific Instruments from Ft. Collins, Colorado, in the 


southern test region that is shown in Figures 3 and 
4. This network is designed for testing the runoff 
and environmental effect of cloud-seeding operations. 
It also provides the ideal starting area for the expan- 
sion of SNOPSET from well to less well instrumented 
regions. Within this newly instrumented area also 
exist some old gages (see Fig. 3). The large ratio 
of planned versus existing gages clearly illustrates 
the large instrumentation requirement that would 
exist if parametric models were to replace the 
present statistical discharge prediction models. 

The Skywater test site consists of a number of 
small drainage basins. One of these, the Wolf 
Creek Drainage Basin, was chosen as the initial 
starting point for the test of SNOPSET. This basin 
is shown in Figure 4. It varies in elevation from 
about 8700 ft to slightly more than 12 000 ft. It 
receives approximately 600 in. of precipitation annu- 
ally, covers approximately 14.5 square miles, and 
is approximately 80 percent forested. It has steep 
slopes facing generally southeast and northwest. It 
is generally snow -covered from mid-October to mid- 
June and has past records of stream gage and pre- 
cipitation and temperature. 

Many instrumentation problems exist, such as 
freezing of ink and storm damage. To avoid such 
loss in input continuity, someone must regularly 
service and maintain all gages. This is a severe 
operations problem in the high mountain areas that 
hold most of the snowpack. To illustrate, from all 
existing old gages shown in Figure 3, only the single 
location near the summit of Wolf Creek Pass was 
kept operating continuously enough to support stream- 
flow prediction. Clearly, the cost of developing and 
maintaining an adequate precipitation network for 
parametric stream discharge prediction in the main 
precipitation area of Figure 2 would be high. We 
hope that this cost might be reduced by augmenting 
a much smaller number of gages with remote obser- 
vations. 

Discharge Prediction Without 
Remote Observations 

The Wolf Creek Drainage Basin has only a single 
location where the continuity of old precipitation 
readings is sufficient for parametric models. How- 
ever, the use of only one precipitation gage is usually 
not sufficient to adjust the model parameters for a 
representation of the drainage basin at hand. The 
adjustment of SNOPSET was nonetheless attempted 
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by augmenting precipitation with temperature records 
(daily maximum and minimum) . 

Figure 5 indicates the results from the first 
guess at the hydrological parameters in the model 
using the 1968-1969 water year as a test case. The 
correlation between recorded and synthesized mean 
daily streamflows is quite poor. Figure 6 shows the 
same results after 19 self-calibrating iterations of 
SNOPSET. The correlation in this case is very good 
(0. 94) . The associated choice of model parameters 
was then verified by applying the unchanged model to 
the 1969-1970 water year. 

The results are shown in Figure 7. It is appar- 
ent that the "plant transfer function" has been ade- 
quately determined by application of SNOPSET. 

The good correlation between computer and pre- 
dicted stream discharge indicates that the use of 
parametric models might offer opportunities for re- 
ducing the density of streamflow gages. However, 
Figure 7 does not imply any forecast ability since 
model input and output overlap in time. 

The predictive ability of SNOPSET was tested by 
extrapolating the streamflow beyond the time of the 
precipitation and temperature inputs. The results 
are shown in Figure 8. The abscissa shows the fore- 
cast period; i. e., the period between the computed 
discharge and the last precipitation measurement. 

The ordinate shows the quality of the prediction in 
terms of a standard deviation between the daily aver- 
ages of predicted and measured streamflow. This 
deviation is approximately 3 times the residual devi- 
ation that remained after adjusting the SNOPSET 
parameters as shown in Figure 7. Apparently, the 
model provides a good prediction over a forecast 
period of approximately 6 months. This is adequate 
for covering the period of snowmelt. During these 
forecasts, the predicted and measured streamflows 
correlate within approximately 75 percent. For a 
forecast in excess of 6 months, the deviation between 
measurement and forecast increases more and more. 
We believe that the above results warrant to base a 
continuation of our present demonstration project and 
suitable updated versions of SNOPSET, that would 
take remote observation in addition to readings from 
a few precipitation and temperature gages. 


Aircraft Instrumentation for 
Augmentation of Precipitation 
Measurements 

Our first demonstration on the interpolation of 
precipitation gage readings with remote observations 
will be restricted to the snowpack in high mountain 
areas. The area extent of snow can easily be detect- 
ed from space and the depletion of the snow cover 
appears to be directly related to runoff characteris- 
tics [ 3 through 10] . Even the low resolution of the 
meteorological satellites gave adequate estimates of 
the yearly precipitation input in the Sierra Nevadas. 
Monthly inputs to smaller drainage areas can hope- 
fully be derived from repetitive surveys that provide 
higher resoultion images, such as the Earth Resources 
Technology Satellite (ERTS-A). Additional significant 
information for discharge prediction can hopefully be 
derived from incremental changes of the melt line. 

The boundary of melting snow should be accessible 
from the 0° C contour of a thermal map. 

Any application of space observations to stream 
discharge prediction would lose most of its value if 
the precipitation estimate does not include rain be- 
sides snow. Unfortunately, the interpolation of rain- 
fall gages is much more difficult than the above ex- 
trapolation of snowfall measurements since soil 
moisture is much more difficult to detect. 

Temperature anomalies and greening of dry vege- 
tation have been proposed as qualitative and indirect 
indicators of rainfall. However, then indicators may 
have many other causes besides recent rain and are 
thus probably not feasible for the interpolation of 
rainfall gages. We propose instead to use a map of 
the vertically integrated water vapor mass. Incre- 
mental changes of such a map should provide infor- 
mation on the water loss by evaporation and transpi- 
ration, which is also indicated by the decrease of the 
precipitation level in the rain gages. Repetitive sur- 
veys of the water vapor mass distribution should thus 
allow interpolation of rain gages shortly after precipi- 
tation inputs when these levels decline. Hopefully, 
this interpolation will also hold for the immediately 
preceding period of rainfall that was indicated by 
rising gage levels. 

The CSU aircraft with a special scanner to pro- 
vide maps of snow cover, temperature, and vertically 
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integrated water vapor mass ( Fig. 9) is being in- 
strumented by MSFC. This scanner has three chan- 
nels, two for reflected sunlight (0. 83 to 0. 87Mm and 
0. 91 to 0. 93*m) , and one for thermal emission (8 to 
12Mm) . The restriction to only three channels was 
necessary to conserve weight and funds so that a 
small aircraft may be used for surveying hourly and 
daily variations of precipitation as required by the 
direct surveillance network. 

The bandpass of the first reflectance channel is 
chosen so that atmospheric propagation effects of 
water vapor are minimized. This channel will be 
used to provide the snow maps. A normalizing fac- 
tor is measured directly per scan line by viewing 
the above sky screen. By using this additional ref- 
erence pulse, channel 1 can be used to correct for 
extraneous illumination factors such as partial cloud 
cover. 

The thermal channel is viewing a sky screen 
once per scan mirror revolution to obtain a signal 
that is related to the ambient temperature of the 
aircraft housing. In addition, the thermal channel 
will record the emission of an adjustable and stabi- 
lized blackbody calibration source. These two ad- 
ditional calibration pulses are then used to convert 
the thermal signal into a radiometric temperature 
that should provide a good approximation of the sur- 
face temperature. 

The bandpass of the second reflectance channel 
coincides with the absorption band of water vapor at 
0. 93p. This channel should thus be used to map the 
transmission loss that is provided by the atmospheric 
water vapor. 

This transmission must then be interpreted in 
terms of integrated water vapor mass (precipitable 
centimeters) . Background for such spectroscopic 
interpretation is available from a 15-year informal 
and international cooperation on the propagation 
characteristics of the water vapor molecule. The 
results of this investigation are summarized by the 
Wave Propagation Laboratories of the National 
Oceanic and Atmospheric Administration (NOAA) 

£ 16] . Their Slant-Path computer codes have been 
made available to MSFC as part of a previous joint 
field test program. 

The proposed water-vapor map would account for 
the total water-vapor mass between the ground and 
the scanner; whereas, the hydrological applications 
refer only to the portion that is related to the evapo- 
ration from the underlying soil and the transpiration 


from plants . The balance is provided by humidity 
in converted air masses. As long as this convection 
circulates the air within the large drainage area of 
interest, one might expect that its effect would can- 
cel out in the integration across this area [ 17] . 
However, a significant influx of atmospheric mois- 
ture into the drainage area might present interpre- 
tation problems. In such an event, statistical cor- 
relation concepts are conceivable which might re- 
trieve the evapotranspiration component by a digital 
correlation of water-vapor and temperature maps. 

Conclusion 

The successful statistical prediction of stream 
discharge from historical records of a few gages im- 
plies that the spatial interpolation of precipitation 
gages does not need high local accuracy for obtaining 
acceptable overall precipitation inputs. Remote ob- 
servations might thus have a chance for spatial inter- 
polation of a few gages over large regions, as shown 
in Figure 2. The above successful forecast over 6 
months also implies that the temporal interpolation 
between successive overflights does not need a great 
local accuracy. Temporal interpolation errors are 
obviously acceptable within the deviations of daily 
averages that are indicated by the statistical hydro- 
logical models in present use. A 9-day interval be- 
tween remote observations of half the drainage area 
is within the capability of existing unmanned satel- 
lites (ERTS) and might suffice for the temporal 
interpolation. 

A special aircraft scanner is being built and 
tested which should provide maps on the aerial ex- 
tent of snow, on surface temperature and on the ver- 
tically integrated water vapor mass for the spatial 
and temporal interpolation of a few precipitation 
gages. Such demonstrations should be conducted in 
large drainage areas where statistical hydrological 
models are being used. The errors of the local in- 
terpolation can be spot checked in small subregions, 
where the existing gage network is sufficient to sup- 
port the adjustment of the SNOPSET parameters for 
local streamflow estimates without remote observa- 
tions. The errors of the local space and time inter- 
polation would thus be established by reducing the 
number of gage inputs while simultaneously using 
remote observations. The errors of the overall 
precipitation inputs and evapotranspiration losses 
would subsequently be established by estimating the 
propagation and self- cancellation of the local inter- 
polation errors in the space and time integrals that 
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establish the overall precipitation inputs to the 
drainage area. The feasibility of using space obser- 
vations for precipitation estimates becomes apparent 
by comparing the estimated overall errors with the 
tolerances that were established when using the 
statistical forecast method. If feasible, the use of 
space observations would then provide for adjusting 
the statistical methods for changes in the drainage 
area by cloud-seeding operations, new river chan- 
nels, new reservoirs, etc. Equally, if not more 
important, the use of space observations might pro- 
vide stream discharge forecasts in undeveloped 
regions, where several decades of precipitation and 
discharge records do not yet exist. 

Evaporation losses usually exceed the water 
flowing in the tributaries and their accurate knowl- 
edge is thus necessary if one wants to predict dis- 
charge by subtracting precipitation and evaporation. 
Most parametric hydrological models estimate evap- 
oration losses with empirical factors which account 
for various soils, vegetation covers, slopes, etc. , 
in the drainage basin. Many of these factors might 
not be needed if the evaporation loss can be estimat- 
ed from the water vapor mass. The survey of at- 
mospheric water vapor distribution is thus not only 
needed for estimating overall rainfall, but may also 
assist in developing parametric stream discharge 
prediction models in sparsely instrumented drainage 
areas. 
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Figure 1. Hydrological applications in Colorado and Tennessee Valley. 



Figure 2. Reservoir operations by Bureau of Reclamation. 
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Figure 7. Verification of hydrological model through 1969-1970 stream discharge prediction. 


69/70 DEVIATION BETWEEN THE FORECAST AND THE 
DAILY DISCHARGE MEASURE AVERAGED OVER THE YEAR 



Figure 8. Streamflow prediction with self-calibrating extension at Kentucky Watershed Model. 
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